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We report a measurement of the production ratio of charged and neutral B mesons from T(4S f ) 
decays based on the ratio of efficiency-corrected yields for the charmonium modes J/ip K + and J/ip K° 
with 81.9 fb -1 of data collected with the BABAR detector on the T(4S) resonance at 10.580 GeV. We 
find a value of 1.006±0.036(s£at) ±0.031(sys) for the ratio R +/0 = T(Y(4S) -> B+ B~) /Y(T{AS) -> 
B°B°). 
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PACS numbers: 13.25.Hw, 13.25.Gv, 14.40.Nd 

A measurement of the B + /B° production ratio 

R+/0 = T(T(AS)^B+B-) 
T(r(4S) B°B°) 

from the T(AS) meson is an essential element in deter- 
mining branching fractions and quark-mixing matrix el- 
ements at the B factory experiments. It can also provide 
information about the structure of the T(4S') meson that 
can be used to discriminate between available models. 

Over the past 15 years it has been frequently assumed 
that is equal to one, although many models predict 

that this may not be the case. Early calculations pre- 
dicted that the ratio could be up to 20% greater than one, 
due to large Coulomb corrections Q . Taking into account 
the structure of the B and T(4S) reduces the effect of the 
Coulomb interaction and can even lead to the ratio being 
less than unity [2J . With the prospect of precision mea- 
surements from the B factories, there has been a recent 
revival in theoretical work on the subject. A more de- 
tailed calculation has been done in a non-relativistic effec- 
tive field theory with B* intermediate states in the pion 
potential, which introduces isospin-breaking in strong in- 
teractions. These calculations predict a value 1.1 — 1.2 0. 
Other calculations attempting to take into account the 
structure of the mesons and hadronic final state interac- 
tions predict a ratio 0.9 — 1.2 4], but with rapid variation 
as a function of the center-of-mass energy near the T(4S) 
resonance. However, such rapid variation in the charged- 
to-neutral ratio has not been seen in scans across the 
0(1020) resonance |5J. For the T(AS), there are published 
measurements of R+/° by CLEO (1.04 ± 0.07 ± 0.04 @, 
1.058±0.084±0.136 0), BABAR (1.10±0.06±0.05 with 
20ffcT 1 ), and Belle (1.01 ± 0.03 ± 0.09 9] with 29fb~ 1 ). 
Now that a significantly larger Y(AS) data sample is 
available at BABAR we can reduce the statistical uncer- 
tainty to the point where it is possible to confront the 
various theoretical predictions. 

In this analysis we use the decay modes B° — ► J/ipK® 
and B+ -> J/ipK + 10], where J/ip -► £+£' and K° s -> 
7r + 7r~, to measure the B + /B° production ratio. These 
decays are good candidates for measuring R + / Q since 
isospin violation in the B — ► J/ip K decays is expected to 
be small in the Standard Model, of order A 3 « 0.01 
when rescattering is small, where A is defined as the sine 
of the Cabibbo angle. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II e+e~ storage ring. The 
data sample corresponds to 81.9 fb _1 of integrated lumi- 
nosity collected at the T(45) resonance. The distribution 
of center-of-mass energies due to the beam-energy spread 
is Gaussian with a = 4.6 MeV [12j . The mean energy 
of our sample is 10.580 GeV, with all data accumulated 
within one sigma of this value. 



The BABAR detector is fully described elsewhere [l3| . 
It consists of a charged-particle tracking system, a 
Cherenkov detector (DIRC) for particle identification, 
an electromagnetic calorimeter, and a system for muon 
identification. The tracking system consists of a 5-layer, 
double-sided silicon vertex tracker and a 40-layer drift 
chamber (filled with a mixture of helium and isobutane), 
both in a 1.5-T magnetic field supplied by a supercon- 
ducting solenoidal magnet. The DIRC is an imaging 
Cherenkov detector relying on total internal reflection in 
the radiator. The electromagnetic calorimeter consists of 
6580 CsI(Tl) crystals. The iron flux return is segmented 
and instrumented with resistive plate chambers for muon 
identification. 

Hadronic events are selected by requiring the presence 
of at least three tracks in the angular region 0.41 < 
$lab < 2.54 rad, where Slab is the polar angle with 
respect to the beam direction. The ratio between the 
2nd and 0th order Fox- Wolfram [14| moments must be 
less than 0.5. We also require that the total energy of 
all particles in the event be greater than 4.5 GeV. The 
primary vertex, which is constructed from charged tracks 
with impact parameter less than 1mm in the plane trans- 
verse to the beam direction, must be within 0.5 cm of the 
beam spot in the plane transverse to the beam direction 
and within 6 cm along the beam direction. 

We reconstruct candidates for J/ip mesons in the decay 
modes J/ip — > e + e _ and For J/ip — > e + e~ decays 

one track is required to pass a tight electron selection and 
the other a loose requirement [8|, while for J/ip — * /i + 
decays we require one track to pass a loose muon selec- 
tion and the other a minimum-ionizing requirement ||. 
The daughter tracks of the J/ip candidate are required 
to have 12 hits in the drift chamber, lie in the angu- 
lar range 0.41 < 6*lab < 2.409 rad for electrons and 
0.41 < 6>lab < 2.54 rad for muons, and have a trans- 
verse momentum of at least 100 MeV/ c. To increase the 
efficiency of the event selection, electron candidate tracks 
are combined with photon candidates to recover some of 
the energy lost in bremsstrahlung || . A geometric vertex 
constraint fit is applied to the lepton track pair. The in- 
variant mass requirements for the J/ip —> e + e~ and J/ip 
— > [i+ir channels are 2.95 < M e + e - < 3.14GeV/c 2 and 
3.06 < Af M+M - < 3.14GeV/c 2 . We require |cos^| to be 
less than 0.8 and 0.9 for J/ip — » e + e~ and J/ip — > 
respectively. The helicity angle 6t is the angle in the J/ip 
rest frame between the positively charged J/ip daughter 
and the reversed K flight direction in the B meson rest 
frame. 

We reconstruct Kg meson candidates from two charged 
tracks, which are not required to originate from the in- 
teraction point or to have drift chamber hits, in contrast 
to the J/ip daughters. The tracks are assigned the pion 
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mass to compute M 7r + 7r - , which is required to lie in the 
range 0.490 — 0.505 GeV/c 2 . Also, in order to reject com- 
binatorial background, we only retain candidates with a 
fitted K® vertex displaced more than 1 mm from the J/0 
vertex. Candidates for K + mesons are assigned the kaon 
mass and are required to form a vertex with the J/0 can- 
didate. No particle identification requirements are made 
for this track. 

The selection of B candidates relies on the kinematic 
constraints given by the T(4S) initial state. Two largely 
uncorrelated variables are used: the energy-substituted 
B mass toes = \J ( s /2 + Po ■ Pb) 2 /-^ 2 ~~ Vbi wnere the 
subscripts and B refer to the e + e~ system and the 
B candidate respectively, s is the square of the center- 
of-mass energy, and energies (E) and momentum vec- 
tors (p) are computed in the laboratory frame; and 
AE = E* B — y/s/2, where Eg is the B candidate energy 
in the center-of-mass frame. In cases where multiple B 
candidates are present in the same event, 2% of the total, 
only the one with the smallest absolute value of AE is 
retained. 

The signal region in the TOes-A-E plane is defined by 
5.27 < toes < 5.29GeV/c 2 and \AE\ < 3a(AE). The ob- 
served resolutions for data and Monte Carlo for the differ- 
ent modes are listed in Table [I] The toes sideband is de- 
fined by 5.20 < toes < 5.27GeV/c 2 and \AE\ < 3a(AE). 
Upper and lower AE sidebands, used for the evaluation 
of systematic errors, are defined as 50 < AE < 120 McV 
and -120 < AE < -50MeV. 

Since we are measuring the ratio of B + to B° 
efficiency-corrected yields many of the selection require- 
ments are in common and have been optimized previ- 
ously |8|. Therefore for this analysis, we have only re- 
considered the optimization of the statistical uncertainty 
of the measurement due to those requirements that are 
different for the two modes. These requirements include 
the K° flight length, K° mass window, and the AE win- 
dow. The optimization of these variables maximizes the 
ratio N cand /y/a 2 (N cand ) + a 2 (N hkg ) where cr(A cand ) and 
c(iVbkg) are the uncertainties on the number of signal 
candidates N cand predicted by Monte Carlo (MC) simu- 
lation and combinatorial background iVbkgj respectively, 
that pass the event selection procedure. iV C and is defined 
as number of events in the signal region. 

TABLE I: Summary of the resolution for AE in data and MC 
simulation. 



Mode a(AE) [MeV] 

B J/0 MC Data 

B + e + e~ 11.42 ±0.11 10.87 ± 0.25 
B+ n+fj- 9.72 ±0.07 9.25 ±0.20 
B u e+e" 9.50 ±0.11 10.02 ± 0.42 
B° n + fi~ 7.92 ±0.07 8.52 ±0.32 

We fit the toes distribution in the toes sideband with 



an empirical phase-space-motivated function introduced 
by ARGUS [l^. The fitted distribution is then inte- 
grated over the signal region to determine the number 
of combinatorial background events A^kg- I n addition to 
combinatorial backgrounds there are other background 
sources, mostly in B decays to charmonium, that peak 
near the B mass in toes- These peaking sources are 
negligible for the neutral B sample, but include small 
contributions from B° — > J/0 K® and B + — > J/07T+ for 
the charged B sample. Requiring particle identifications 
on the K + candidate will reduce these contributions but 
introduces a larger systematic error. To determine the 
number of background events that peak in the toes signal 
region, A pea k, we use appropriately combined MC sam- 
ples of continuum e + e~ and generic [3 BB events (with 
signal events removed), which have been scaled to the 
integrated luminosity of the data sample. This distribu- 
tion is then fitted with an ARGUS function as described 
above. We determine A peak by counting the number of 
events in the signal region and subtracting the integral of 
the ARGUS function over this same region. The signal 
yield is then defined by A signa i = A cand - A bkg - A peak . 
The observed distributions in toes an d AE for B° — > 
J/ipK® and B + — > J/ipK + candidates in data are shown 
in Figs. Q] an d|21 respectively. 

The efficiency-corrected ratio of observed events is 
given by: 

ggnal/fj = +/0 g(i? + ^ J/0g+) 
^s°i g n al /(Ao) B{B*^ JftK*) 

= +/0 2T(B+^J/^K+)r + 
T(B° -> J/iPK°)t 

where / = 68.60 ±0.27% I7!j is the K2 -> tt+tt" branch- 
ing fraction, t+/t = 1.083 ±0.017 [Hf is the ratio of B+ 
and £>° lifetimes, and e is the selection efficiency. There- 
fore, assuming isospin invariance in the B — > J/ipK de- 
cay, T(B+ -> J/0A'+) = T(B° -> J/tpK°) 18], the ratio 
of efficiency-corrected yields is determined from: 

+ /0 jVsignalfo/ TO , , 

R (1) 

The ratio of efficiency-corrected yields is determined sep- 
arately for J/0 — > e + e~ and J/0 — > so that lepton 
identification efficiencies cancel. The separate measure- 
ments are then averaged, keeping track of correlated un- 
certainties, to produce a final value for R + /° . 

Sources of systematic uncertainties can be classified 
into those arising from uncertainties on efficiencies and 
those from candidate selection and backgrounds. The 
efficiency uncertainties are due to K® reconstruction, 
tracking, and kaon/pion tracking efficiency differences. 
In the ratio of the efficiency-corrected yields, the track- 
ing uncertainty is due to the extra track required to re- 
construct the B° — > J/ipK® mode. We determine the 
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FIG. 1: Distribution of m ES for \AE\ < 3cr in the B° -> 
J/0 Ks sample for (a) J/0 — > e + e~ and (b) J/tp — > /Lt + H~ ■ 
Distribution of AE for m E s > 5.27GeV/c 2 in the 5° -> 
J/ipKg sample for (c) J/ip — > e + e~ and (d) J/0 — > ^ + ^ ■ 



FIG. 2: Distribution of m E s for |A£| < 3a in the B + -> 
J/ip K + sample for (a) J/0 — ► e + e~ and (b) J/-!/ 1 — ► 
Distribution of AS for m E s > 5.27GeV/c 2 in the B+ -> 
J/i/)AT+ for (c) J/0 -> e+e~ and (d) J/0 -> n + n~- 



relative kaon/pion tracking reconstruction efficiency by 
comparing the ratio of efficiencies for B + — ► J/0 AT + and 
_B + — » J/0 7r + Monte Carlo. The systematic error of 
0.6% is taken to be half the size of the estimated differ- 
ence. Finally, for the uncertainty on the K® efficiency 
we take a sample of inclusive K® candidates that are 
binned in transverse momentum (pt), laboratory polar 
angle (Slab), an d transverse flight length (dr). A relative 
correction for reconstruction of a displaced K® candidate 
is determined in each pt and Slab bin by assuming the 
tracking efficiency for a short-lived K® close to the inter- 
action region is the same as for prompt tracks. Thus, the 
ratio of data to MC relative efficiency is normalized to 
unity for small dr and then used to derive a MC correc- 
tion factor for larger displacements. By varying the size 
of the dr, pr, and Slab bins we determine a systematic 
uncertainty for this procedure. The normalization bin 
for the correction is well inside the radius of the beam 



pipe. We vary the definition of this bin as a check of the 
hypothesis that these tracks have the same efficiency as 
normal charged tracks. 

TABLE II: Summary of the relative systematic uncertainties 
on the efficiency-corrected yields. 



Mode Parameters (%) 

ARGUS Peaking Va^ 

B J/0 £Trk /tt+ e K°, Bkgd. Bkgd. Selection Total 

B+ e + e- - 06 05 O O OlT 

B+ ij+jj,- - 0.6 - 0.6 0.4 1.0 1.4 

B u e+e" 1.3 - O 08 O (U 2A~ 

B° n+n' 1.3 - 1.8 0.5 0.1 1.3 2.6 

The selection and background systematic uncertainties 
are attributed to the selection criteria, the ARGUS back- 
ground shape, and the peaking background subtraction. 
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The selection requirements on the K® mass, K® flight dis- 
tance, and AE are varied within reasonable ranges. The 
K° mass and AE selection windows were increased and 
decreased from the nominal value by half a sigma and the 
K® vertex displacement requirement was removed. The 
largest difference from the nominal efficiency-corrected 
yield is taken conservatively as a systematic uncertainty. 
The continuum background systematic uncertainty is de- 
termined by varying the ARGUS parameter by two sigma 
to account for any model dependence. The peaking back- 
ground uncertainty is evaluated from the discrepancy be- 
tween data and MC in the the upper and lower AE side- 
bands. The larger of the two discrepancies is taken as 
the systematic uncertainty. This is a more conservative 
approach than using the uncertainties for the relevant 
branching fractions. Table ITU summarizes the sources of 
systematic uncertainty for this analysis. 



TABLE III: Summary of values needed to determine the effi- 
ciency corrected yields. 



Mode 




Parameters 




B J/tj} 


■^V cand 




Efficiency (%) 


B+ e+e" 


2213 


19.5 ± 5.0 9.6 ± 3.2 


40.8 ± 0.4 


B° e+e- 


502 


2.6 ± 2.0 2.4 ± 1.5 


29.9 ± 0.4 


B+ 


2497 


50.6 ± 7.2 33.5 ± 4.6 


47.8 ± 0.4 


B° n+fT 


577 


2.0 ± 1.5 2.4 ± 2.1 


35.6 ± 0.4 



are not restricted to a small range in \AE\. It also em- 
ploys the MC predicted background distributions, rather 
than the empirical shape imposed by the ARGUS func- 
tion. The non-parametric fit method finds results that 
are consistent with the simpler counting method, both 
for the full sample and for data subsets. 

The observed value for R + /° is close to one, as has 
been assumed by most branching fraction measurements 
obtained on the T(AS), with a ratio as large as 1.2 dis- 
favored at the four sigma level. Our measurement will 
aid in restricting models of T(45) decays. It also al- 
lows a quantitative determination of the contribution 
from to all branching fractions that are deter- 
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NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
INFN (Italy), FOM (The Netherlands), NFR (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
viduals have received support from the A. P. Sloan Foun- 
dation, Research Corporation, and Alexander von Hum- 
boldt Foundation. 



Table liTTl lists the efficiencies, background composition, 
and number of events in the signal region based on the 
one-dimensional fit with a 3(er) AE requirement. Based 
on Eq.(^) we determine: 

R +/0 (e + e~) = 1.019 ± 0.054(stai) ± 0.031 (sys) 
i? +/0 (/i + /O = 0.994 ± 0M9(stat) ± 0.033(sys) 
R +/a (avg) = 1.006 ±0.036(stot)±0.031(sys) 

when assuming isospin conservation in B — > J/tpK de- 
cays. The data sample has a mean energy of 10.580 GeV 
and does not have sufficient spread to test the hypothesis 
of an energy dependent production ratio. 

We have confirmed that the result for the individual 
efficiency-corrected signal yields for the J/ip — > e + e~ and 
J/ip — > M + channels is consistent among seven equal 
subsets of the full sample, as is the ratio of e + e~//i + /i~. 

To check our fitting technique we have performed a 
two-dimensional non-parametric fit to the data. This is 
done by fitting the data to a sum of contributions from 
five different sources (e + e~ — > qq, e + e~ — > cc, generic 
T(4S) -> £ S°generic T(4S*) -> B+B~, and signal) 
whose densities |l9j in AE and toes are determined from 
a non-parametric fit to candidates from Monte Carlo 
samples. The two-dimensional fit is done in the region 
5.200 < toes < 5.270 GeV/c 2 and 0.030 < \AE\ < 
0.120 GeV. This technique has the advantage that we 
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